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Abstract 

The dwarf spheroidal satellite galaxies of the Milky Way are some of the most dark-matter- 
dominated objects known. Due to their proximity, high dark matter content, and lack of as- 
trophysical backgrounds, dwarf spheroidal galaxies are widely considered to be among the most 
promising targets for the indirect detection of dark matter via 7 rays. Here we report on 7 -ray 
observations of 25 Milky Way dwarf spheroidal satellite galaxies based on 4 years of Fermi Large 
Area Telescope (LAT) data. None of the dwarf galaxies are significantly detected in 7 rays, and 
we present 7 -ray flux upper limits between 500 MeV and 500 GeV. We determine the dark matter 
content of 18 dwarf spheroidal galaxies from stellar kinematic data and combine LAT observations 
of 15 dwarf galaxies to constrain the dark matter annihilation cross section. We set some of the 
tightest constraints to date on the annihilation of dark matter particles with masses between 2 GeV 
and lOTeV into prototypical standard model channels. We find these results to be robust against 
systematic uncertainties in the LAT instrument performance, diffuse 7 -ray background modeling, 
and assumed dark matter density profile. 

PACS numbers: 95.35.+d, 95.85.Pw, 98.52.Wz 
Keywords: dark matter; gamma rays; dwarf galaxies 
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I. INTRODUCTION 


Astrophysical evidence and theoretical arguments suggest that non-baryonic cold dark 
matter constitutes ~ 27% of the contemporary Universe [1]. While very little is known about 
dark matter beyond its gravitational influence, a popular candidate is a weakly interacting 
massive particle (WIMP) [2-4], From an initial equilibrium state in the hot, dense early 
Universe, WIMPs can freeze out with a relic abundance sufficient to constitute much, if not 
all, of the dark matter. Such WIMPs are expected to have a mass in the GeV to TeV range 
and a s-wave annihilation cross section usually quoted at ~3 x 1CU 26 cm 3 s -1 (though, as 
pointed out by Steigman et al. [5], a more accurate calculation yields ~5 x lCU 26 cm 3 s -1 at 
masses < 10 GeV and ~2 x 10” 26 cm 3 s -1 at higher masses). In regions of high dark matter 
density, WIMPs may continue to annihilate into standard model particles through processes 
similar to those that set their relic abundance. 

Gamma rays produced from WIMP annihilation, either mono-energetically (from direct 
annihilation) or with a continuum of energies (through annihilation into intermediate states) 
may be detectable by the Large Area Telescope (LAT) on board the Fermi Gamma-ray 
Space Telescope {Fermi) [6]. These 7 rays would be produced preferentially in regions of 
high dark matter density. The LAT has enabled deep searches for dark matter annihila- 
tion in the Galactic halo through the study of both continuum and mono-energetic 7-ray 
emission [7-13]. Additionally, robust constraints have been derived from the isotropic 7-ray 
background [14, 15], galaxy clusters [16], Galactic dark matter substructures [17-19], and 
the study of nearby dwarf spheroidal galaxies [20-24], 

The dwarf spheroidal satellite galaxies of the Milky Way are especially promising targets 
for the indirect detection of dark matter annihilation due to their large dark matter content, 
low diffuse Galactic 7-ray foregrounds, and lack of conventional astrophysical 7-ray produc- 
tion mechanisms [25, 26]. An early analysis by Abdo et al. [20] set constraints on the dark 
matter annihilation cross section from 11-month LAT observations of 8 individual dwarf 
spheroidal galaxies. In a subsequent effort, Ackermann et al. [21] presented an improved 
analysis using 2 years of LAT data, 10 dwarf spheroidal galaxies, and an expanded statis- 
tical framework. The analysis of Ackermann et al. [21] has two significant advantages over 
previous analyses. First, statistical uncertainties in the dark matter content of dwarf galax- 
ies are incorporated as nuisance parameters when fitting for the dark matter annihilation 
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cross section. Second, observations of the 10 dwarf galaxies are combined into a single joint 
likelihood analysis for improved sensitivity. The analysis of Ackermann et al. [21] probes the 
canonical thermal relic cross section for dark matter masses < 30 GeV annihilating through 
the bb and t + t~ channels, while similar studies by Geringer-Sameth and Koushiappas [22] 
and Mazziotta et al. [23] yield comparable results. 

The present analysis expands and improves upon the analysis of Ackermann et al. [21]. 
Specifically, we use a 4- year y-ray data sample with an extended energy range from 500 MeV 
to 500 GeV and improved instrumental calibrations. We constrain the y-ray flux from all 25 
known Milky Way dwarf spheroidal satellite galaxies in a manner that is nearly independent 
of the assumed y-ray signal spectrum (Section III). We use a novel technique to determine 
the dark matter content of 18 dwarf spheroidal galaxies by deriving prior probabilities for the 
dark matter distribution from the population of Local Group dwarf galaxies (Section IV). 
This allows us to increase the number of spatially-independent dwarf galaxies in the joint 
likelihood analysis to 15 (Section V). We develop a more advanced statistical framework to 
examine the expected sensitivity of our search. Additionally, we model the spatial y-ray 
intensity profiles of the dwarf galaxies in a manner that is consistent with their derived 
dark matter distributions. We perform an extensive study of systematic effects arising from 
uncertainties in the instrument performance, diffuse background modeling, and dark matter 
distribution (Section VI). As a final check, we perform a Bayesian analysis that builds on 
the work of Mazziotta et al. [23] and yields comparable results (Section VII). No significant 
y-ray emission is found to be coincident with any of the 25 dwarf spheroidal galaxies. Our 
combined analysis of 15 dwarf spheroidal galaxies yields no significant detection of dark 
matter annihilation for particles in the mass range from 2 GeV to 10 TeV annihilating to 
e + e _ , t + t~, uu, bb, or W + W~ (when kinematically allowed), and we set robust upper 

limits on the dark matter annihilation cross section. 


II. DATA SELECTION AND PREPARATION 

Since the start of science operations in August of 2008, the LAT has continuously scanned 
the y-ray sky in the energy range from 20 MeV to > 300 GeV [6] . The LAT has unprece- 
dented angular resolution and sensitivity in this energy range, making it an excellent instru- 
ment for the discovery of new y-ray sources [27]. We select a data sample corresponding 
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to events collected during the first four years of LAT operation (2008-08-04 to 2012-08-04). 
We use the P7REP data set, which utilizes the Pass-7 event reconstruction and classifica- 
tion scheme [28], but was reprocessed with improved calibrations for the light yield and 
asymmetry in the calorimeter crystals [29]. The new calorimeter calibrations improve the 
in-flight point-spread function (PSF) above ~3GeV and correct for the small (~ 1% per 
year), expected degradation in the light yield of the calorimeter crystals measured in flight 
data. Consequently, the absolute energy scale has shifted upward by a few percent in an 
energy- and time-dependent manner. In addition, the re-calibration of the calorimeter light 
asymmetry leads to a statistical re-shuffling of the events classified as photons. 

We select events from the P7REP CLEAN class in the energy range from 500 MeV to 500 GeV 
and within a 10° radius of 25 dwarf spheroidal satellite galaxies (Figure 1). The CLEAN event 
class was chosen to minimize particle backgrounds while preserving effective area. At high 
Galactic latitudes in the energy range from 1 GeV to 500 GeV, the particle background con- 
tamination in the CLEAN class is ~ 30% of the extragalactic diffuse y-ray background [28] , 
while between 500 MeV and 1 GeV the particle background is comparable to systematic 
uncertainties in the diffuse Galactic y-ray emission. Studies of the extragalactic y-ray back- 
ground at energies greater than 500 GeV suggest that at these energies the fractional residual 
particle background is greater than at lower energies [30]. To reduce y-ray contamination 
from the bright limb of the Earth, we reject events with zenith angles larger than 100° and 
events collected during time periods when the magnitude of rocking angle of the LAT was 
greater than 52°. 

We create 14° x 14° regions-of-interest (ROIs) by binning the LAT data surrounding each 
of the 25 dwarf galaxies into 0.1° pixels and into 24 logarithmically-spaced bins of energy 
from 500 MeV to 500 GeV. We model the diffuse background with a structured Galactic 
y-ray emission model ( gll_iem_v05.fit ) and an isotropic contribution from extragalactic y 
rays and charged particle contamination (iso_clean_v05.txt) } We build a model of point-like 
y-ray background sources within 15° of each dwarf galaxy beginning with the second LAT 
source catalog (2FGL) [27]. We then follow a procedure similar to that of the 2FGL to 
find additional candidate point-like background sources by creating a residual test statistic 
map with pointlike [27]. No new sources are found within 1° of any dwarf galaxy and the 
additional candidate sources have a negligible impact on our dwarf galaxy search. We use the 


1 http : / /f ermi . gsf c . nasa . gov/ssc/ data/ access/lat/BackgroundModels . html 


P7REP_CLEAN_V15 instrument response functions (IRFs) corresponding to the LAT data set 
selected above. When performing the Bayesian analysis in Section VII, we utilize the same 
LAT data set but follow different data preparation and background modeling procedures, 
which are described in that section. 


III. MAXIMUM LIKELIHOOD ANALYSIS 

Limited 7-ray statistics and the strong dependence of the LAT performance on event 
energy and incident direction motivate the use of a maximum likelihood-based analysis to 
optimize the sensitivity to faint 7-ray sources. We define the standard LAT binned Poisson 
likelihood, 

\ nk p—^k 

c( M ,e |®) = n-h-i-’ (!) 

k 

as a function of the photon data, T > , a set of signal parameters, /x, and a set of nuisance 
parameters, 0. The number of observed counts in each energy and spatial bin, indexed by 
k, depends on the data, n = nfc(D), while the model-predicted counts depend on the input 
parameters, A*, = Afc(/x,0). This likelihood function encapsulates information about the 
observed counts, instrument performance, exposure, and background fluxes. However, this 
likelihood function is formed “globally” (i.e., by tying source spectra across all energy bins 
simultaneously) and is thus necessarily dependent on the spectral model assumed for the 
source of interest. To mitigate this spectral dependence, it is common to independently fit 
a spectral model in each energy bin, j (i.e., to create a spectral energy distribution for a 
source) [31]. This expands the global parameters /x and 6 into sets of independent parameters 
{/A j" and ^ ^ ^ . Likewise, the likelihood function in Equation (1) can be reformulated as a 
“bin-by-bin” likelihood function, 

A W. {«;} " II I v i) • (2) 

3 

In Equation (2) the terms in the product are independent binned Poisson likelihood func- 
tions, akin to Equation (1) but with the index k only running over spatial bins. 

By analyzing each energy bin separately, we remove the requirement of selecting a global 
spectrum to span the entire energy range at the expense of introducing additional parameters 
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into the fit. If the source of interest is bright, it is possible to profile over the spectral 
parameters of the background sources as nuisance parameters in each energy bin. However, 
for faint or undetected sources it is necessary to devise a slightly modified likelihood scheme 
where the nuisance parameters in all energy bins are set by the global maximum likelihood 
estimate, 2 

9 = arg max £(/^, 9 \ V ) . (3) 

e 

Fixing the normalizations of the background sources at their globally fit values avoids nu- 
merical instabilities resulting from the fine binning in energy and the degeneracy of the 
diffuse background components at high latitude. Thus, for the analysis of dwarf galaxies 
our bin-by-bin likelihood becomes 


£({m,}. e I v) = n I ®j) • (4) 

j 

The bin-by-bin likelihood is powerful because it makes no assumption about the global signal 
spectrum. 3 However, it is easy to recreate a global likelihood function to test a given signal 
spectrum by tying the signal parameters across the energy bins, 

c^e\v) = l[jc(^),e j (e)\v j ). (5) 

j 

As a consequence, computing a single bin-by-bin likelihood function allows us to subse- 
quently test many spectral models rapidly. In practice, the signal parameters derived from 
the global likelihood function in Equation (5) are found to be equivalent to those derived 
from Equation (1) so long as couplings between the signal and nuisance parameters are 

small. 

For each of the 25 dwarf galaxies listed in Table I, we construct a bin-by-bin likelihood 
function. Normalizations of the background sources (both point-like and diffuse) are fixed 
from a global fit over all energy bins. 4 The signal spectrum within each bin is modeled by a 
power law ( dN/dE oc E -2 ), and the signal normalization is allowed to vary independently 
between bins. In each bin, we set a 95% confidence level (CL) upper limit on the energy flux 

2 For brevity we assume that all nuisance parameters are globally fit, but it is straightforward to generalize 
to the situation where some nuisance parameters are defined independently in each bin. 

3 The second order spectral dependence arising from the choice of a spectral model within each bin is found 
to be ~2% for our energy bin size. 

4 Including the putative dwarf galaxy source in the global fit has a negligible (< 1%) impact on the fitted 

background normalizations. 10 



from the putative signal source by fixing the parameters of all other bins and finding the 
value of the energy flux where the log-likelihood has decreased by 2.71/2 from its maximum 
(the “delta-log-likclihood technique”) [32, 33]. Figure 2 illustrates the detailed bin- by-bin 
LAT likelihood result for the Draco dwarf spheroidal galaxy. 

In Figure 3, we compare the observed bin-by-bin flux upper limits for all 25 dwarf 
spheroidal galaxies to those derived for 2000 realistic background-only simulations. These 
simulations are generated with the LAT simulation tool gtobssim using the in-flight point- 
ing history and instrument response of the LAT and including both diffuse and point-like 
background sources. We find that the observed limits are consistent with simulations of the 
null hypothesis. We emphasize that these bin-by-bin limits are useful because they make no 
assumptions about the annihilation channel or mass of the dark matter particle. 

While the bin-by-bin likelihood function is essentially independent of spectral assump- 
tions, it does depend on the spatial model of the target source. The dark matter distributions 
of some dwarf galaxies may be spatially resolvable by the LAT [34], Thus, for dwarf galaxies 
for which stellar kinematic data sets exist, we model the y-ray intensity from the line-of- 
sight integral through the best-fit dark matter distribution as derived in Section IV B. Dwarf 
galaxies that lack stellar kinematic data sets are modeled as point-like y-ray sources. Possi- 
ble systematic uncertainty associated with the analytic form and spatial extent of the dark 
matter profile are discussed in Section IV B and Section VI. 


IV. DARK MATTER ANNIHILATION 


The integrated y-ray signal flux at the LAT, (p s ( ph cm 2 s 1 ), expected from dark matter 
annihilation in a density distribution, p(r), is given by 


4>s( AD) 


i M f E ~ . f 

4lr 2m DM Je „ dg 7 7 -La 

<&PP 


P 2 ( r 


1.0. S. 


J-factor 



( 6 ) 


Here, the <Lpp term is strictly dependent on the particle physics properties — i.e., the 
thermally-averaged annihilation cross section, (av), the particle mass, ttidm, and the differ- 
ential y-ray yield per annihilation, dV 7 /df7 7 , integrated over the experimental energy range. 5 

5 Strictly speaking, the differential yield per annihilation in Equation (6) is a sum of differential yields into 
specific final states: AN 1 /dE 1 = J2j-Bf dAy/dE 7 , where Bf is the branching fraction into final state /. 
In the following, we make use of Equation (6) in the context of single final states only. 
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The J-factor is the line-of-sight integral through the dark matter distribution integrated over 
a solid angle, Afh Qualitatively, the J-factor encapsulates the spatial distribution of the dark 
matter signal, while $pp sets its spectral character. Both <f>pp and the J-factor contribute a 
normalization factor to the signal flux. 


A. Particle Physics Factor 

We generate 7-ray spectra for dark matter annihilation with the DMFIT package [35] as 
implemented in the LAT ScienceTools. We have upgraded the spectra in DMFIT from 
Pythia 6.4 [36] to Pythia 8.165 [37]. We generate 7-ray spectra for the e + e~ and ii + /i ~ 
channels from analytic formulae, i.e., Equation (4) in Essig et al. [38] (see also [39, 40]). For 
the muon channel, we add the contribution from the radiative decay of the muon (/i - — > 
e~ P e 17 7) using Equations (59), (61), and (63) from Mardon et al. [41]; this can have a 
factor-of-two effect on the integrated photon yield for mp m ~ 2 GeV. We also extend the 
range of dark matter masses in DMFIT to lower (< 10 GeV) and higher (> 5 TeV) masses [35]. 
Thus, the upgraded DMFIT includes 7-ray spectra for dark matter annihilation through 12 
standard model channels (e + e _ , g + /a~ , t + t~ , uu, dd, cc, ss, bb, tt, W + W~, Z°Z°, gg ) for 24 
dark matter masses between 2 GeV and 10 TeV (when kinematically allowed). Annihilation 
spectra for arbitrary dark matter masses within this range are generated by interpolation. 
Appropriate flags are set in Pythia to generate spectra for dark matter masses below 10 GeV. 
However, we note that we ignore the effect of hadronic resonances, which could be important 
for some dark matter models with masses less than 10 GeV. In addition, the hadronization 
model used in Pythia becomes more uncertain at lower energies. The new spectra are 
consistent with those calculated in [35] for dark matter masses between 10 GeV and 5 TeV. 
Outside of this range we believe that the new spectra are a more accurate representation 
of the underlying physics since they are derived by modeling some of the relevant physical 
processes, rather than by ad hoc extrapolations. 


B. J-Factors for Dwarf Spheroidal Galaxies 

The dark matter content of dwarf spheroidal galaxies can be determined through dynam- 
ical modeling of their stellar density and velocity dispersion profiles (see a recent review 
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by Battaglia et al. [42]). Recent studies have shown that an accurate estimate of the dy- 
namical mass of a dwarf spheroidal galaxy can be derived from measurements of the average 
stellar velocity dispersion and half-light radius alone [43, 44], The total mass within the 
approximate half-light radius and likewise the integrated J-factor have been found to be 
fairly insensitive to the assumed dark matter density profile [45, 46]. We use the prescrip- 
tion of Martinez [47] to calculate J-factors for a subset of 18 Milky Way dwarf spheroidal 
galaxies possessing dynamical mass estimates constrained by stellar kinematic data sets (Ta- 
ble I) [48-55], 

We construct a likelihood function for each individual dwarf spheroidal galaxy from the 
observed luminosity, half-light radius, r^, and mass within the half-light radius, Mh, together 
with their associated uncertainties. We derive priors on the relationships between the lumi- 
nosity, L, maximum circular velocity, R max , and radius of maximum circular velocity, i?y max , 
from the ensemble of Local Group dwarf galaxies [56] . By implementing this analysis as a 
two-level Bayesian hierarchical model [57], we are able to simultaneously constrain both the 
priors and the resulting likelihoods for the dynamical properties of the dwarf galaxies [47]. 
This approach benefits from utilizing all available knowledge of the dark matter distribu- 
tion in these objects while mitigating any systematic differences between the observed dwarf 
galaxies and numerical simulations [58, 59]. 

While priors on the relationships between, L, V max , and i7y max can be derived from the 
data, their form must be chosen a priori. Motivated by simulations [58, 59], we assume 
a prior on the relationship between i?y max and Vmax that follows a linear form between 
log(i?y max ) and log(Wiax) with an intrinsic Gaussian scatter A f(n,cr), 

'P (log(7?y max ) | log(R max )) « N(a rv log(l/ max ) + /3 rv , a rv ). (7) 

Similarly, we assume that the prior in log(Wiax) is linearly related to log(L) with Gaussian 
scatter, 

P(log(V max ) | log(L)) « JV(a v i log (L) + 0 vh a vt ). (8) 

Finally, we assume a power-law prior on the galaxy luminosity function [60], 

V(L)^L aL . (9) 
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The 7 free parameters of these priors, {a rv , f3 rv , a rv , a v i, /3 v i , a v i , ql}, are derived from the full 
Local Group dwarf galaxy sample. Additionally, each dwarf is described by its measured half- 
light radius, the mass within this radius, the total luminosity, and the associated errors on 
these quantities. We constrain the full set of parameters with a Metropolis nested sampling 
algorithm with approximately 500,000 points per chain [61]. 


The inner slope of the dark matter density profile in dwarf galaxies remains a topic of 
debate [62-65]. Thus, we repeat the above procedure for two different choices of the dark 
matter profile: (i) a cuspy Navarro-Frenk- White (NFW) profile [66], and (ii) a cored Burkert 
profile [67]: 


P(r) 

p(r) 


Por, 


r{r s + r) 2 

Po r s 


r s + r)(r 2 + r 2 


NFW 

(10) 

Burkert 

(11) 


After deriving the halo parameters for each form of the dark matter density profile, we 
perform the line-of-sight integral in Equation (6) to calculate the intensity profile and in- 
tegrated J-faetor within an angular radius of 0.5° (Af2 ~ 2.4 x 10~ 4 sr). Values of the 
integrated J-factor and spatial extension parameter (defined as the angular size of the scale 
radius) are included in Table I. We cross check these results with those derived from an anal- 
ysis using priors derived from numerical simulations [21, 45, 68] and flat, “non-informative” 
priors [34, 56]. We find that the resulting J-factors are consistent within their stated uncer- 
tainties. We confirm that the integrated J-factor within 0.5° is fairly insensitive to the choice 
of dark matter density profile so long as the central value of the slope is less than 1.2 [68]. 
Thus, we default to modeling the dark matter distributions within the dwarf galaxies with 
a NFW profile and examine the impact of this choice in more detail in Section VI. 


14 



V. CONSTRAINTS ON DARK MATTER 


A. Individual Dwarf Galaxies 


For each of the 25 dwarf galaxies listed in Table I, we use Equation (5) to create global 
likelihood functions for 6 representative dark matter annihilation channels (e + e _ , /x + /W, 
t + t~, uu, bb, and W + W ~ ) and a range of particle masses from 2 GeV to lOTeV (when 
kinematically allowed). We define a test statistic (TS) as the difference in the global log- 
likelihood between the null (ji = fi 0 ) and alternative hypotheses (fi = fx) [27, 69], 


TS 


-2 In 


£(Mo,8|P) \ 

c(ii,e\v)) 


( 12 ) 


We note that this definition of TS differs slightly from that commonly used, since the back- 
ground parameters are fixed at their best-fit values (Section III). We find that the TS dis- 
tribution derived from simulations is well matched to a y 2 /2 distribution with one bounded 
degree of freedom, as predicted by the asymptotic theorem of Chernoff [70] (Figure 4). How- 
ever, the study of random high-latitude blank fields suggests that the significance calculated 
from asymptotic theorems is an overestimate of the true signal significance (Section VI). 
We find no significant y-ray excess coincident with any of the 25 dwarf galaxies for any 
annihilation channel or mass. 6 

The largest deviation from the null hypothesis, TS ~ 7, occurs when fitting the Sagittar- 
ius dwarf spheroidal galaxy. Sagittarius is located in a region of intense Galactic foreground 
emission and is coincident with proposed jet-like structures associated with the Fermi Bub- 
bles [71, 72], We thus find that systematic changes to the model of the diffuse y-ray emission 
(Section VI) can lead to TS changes of ~ 5 in this region. 

To derive constraints on the dark matter annihilation cross section, we utilize the kinemat- 
ically determined J-factors derived in Section IV B. We select a subset of 18 dwarf galaxies, 
excluding the 7 galaxies that lack kinematically determined J-factors: Bootes II, Bootes III, 
Canis Major, Leo V, Pisces II, Sagittarius, 7 and Segue 2. Statistical uncertainty in the 
J-factor determination is incorporated as a nuisance parameter in the maximum likelihood 

6 We note that J-factors are not required to test the consistency of the y-ray data with spectral models of 
dark matter annihilation. 

7 While some authors have suggested a large dark matter component in the Sagittarius dwarf [73], tidal 
stripping of this system leads to complicated stellar kinematics [74]. 


15 



formulation. Thus, the likelihood function for each dwarf galaxy, i, is described by, 


a i I A) =A(Ab Oi I A) 

X ^ e -( 1 °glo( J i)- 1 °glo( J i)) 2 / 2o ' l 2 _ ^ ^ 

ln( 10 ) Ji\f 7 knOi 

Here, denotes the individual LAT likelihood function for a single ROI from Equation (5) 
and the oli include both the flux normalizations of background 7 -ray sources (diffuse and 
point-like) and the associated dwarf galaxy J-factors and statistical uncertainties. We fix 
the background normalizations at their best-fit values and profile over the J-factor uncer- 
tainties as nuisance parameters to construct a 1 -dimensional likelihood function for the dark 
matter annihilation cross section [33]. We note here that the J-factor uncertainties must be 
incorporated as nuisance parameters at the level of the global maximum likelihood fit and 
not on a bin-by-bin basis. Applying J-factor uncertainties to each energy bin individually 
would multi-count these uncertainties when creating the global likelihood function. 

For each of the 18 dwarf galaxies with kinematically determined J-factors, we create 
global likelihood functions incorporating uncertainties on the J-factors and derive 95% CL 
upper limits on the dark matter annihilation cross section using the delta-log-likclihood 
technique (Tables 11- VII). From simulations of the regions surrounding each dwarf galaxy, 
we confirm the well-documented coverage behavior of this technique, finding that the limit 
often tends to be somewhat conservative [33]. As can be seen from Tables IV- VI, the results 
for several individual dwarf galaxies (i.e., Draco, Coma Berenices, Ursa Major II, Segue 1, 
and Ursa Minor) probe the canonical thermal relic cross section of 3 x 10 - 26 cm 3 s -1 for low 
dark matter masses. 

B. Combined Analysis 

Under the assumption that the characteristics of the dark matter particle are shared 
across the dwarf galaxies, the sensitivity to weak signals can be increased through a combined 
analysis. We create a joint likelihood function from the product of the individual likelihood 
functions for each dwarf galaxy, 


£(/U {«*} I T>) = U A(/U on | Vi) . 

i 


( 14 ) 
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The joint likelihood function ties the dark matter particle characteristics (i.e., cross section, 
mass, branching ratio, and thus 7 -ray spectrum) across the individual dwarf galaxies. From 
the starting set of 18 dwarf galaxies with kinematically determined J-factors, we select 
a set of non-overlapping dwarfs to ensure statistical independence between observations. 
When the ROIs of multiple dwarf galaxies overlap, we retain only the dwarf galaxy with 
the largest J-factor. This selection excludes 3 dwarf galaxies: Canes Venatici I in favor of 
Canes Venatici II, Leo I in favor of Segue 1, and Ursa Major I in favor of Willman 1. Thus, 
we are left with a subset of 15 dwarf spheroidal galaxies as input to the joint likelihood 
analysis: Bootes I, Canes Venatici II, Carina, Coma Berenices, Draco, Fornax, Hercules, 
Leo II, Leo IV, Sculptor, Segue 1, Sextans, Ursa Major II, Ursa Minor, and Willman 1. 

No significant signal is found for any of the spectra tested in the combined analysis, 
prompting us to calculate 95% CL constraints on (av) using the delta-log-likelihood pro- 
cedure described above (Figure 5). We repeat the combined analysis on a set of random 
high-latitude blank fields in the LAT data to calculate the expected sensitivity (Section VI). 
When calculating the expected sensitivity, we randomize the nominal J-factors in accord 
with their measurement uncertainties to form an unconditional ensemble [75, 76]. Thus, 
we note that the positions and widths of the expected sensitivity bands in Figure 5 reflect 
the range of statistical fluctuations expected both from the LAT data and from the stellar 
kinematics of the dwarf galaxies. The largest deviation from the null hypothesis occurs for 
dark matter masses between 10 GeV and 25 GeV annihilating through the bb channel and 
has TS = 8.7. The conventional threshold for LAT source discovery is TS > 25 [27]. Studies 
of random blank-sky locations (Section VI) show that this TS cannot be naively converted 
to a significance using asymptotic theorems. Thus, we derive a global p-value directly from 
random blank high-Galactic-latitude sky positions yielding p « 0.08. This global p-value 
includes the correlated trials factor resulting from testing ~ 70 dark matter spectral mod- 
els. As noted by Hooper and Linden [7], it is difficult to distinguish the 7 -ray spectra of a 
~ 25 GeV dark matter particle annihilating into bb and a lower-mass dark matter particle 
(~5GeV) annihilating to t + t~ . However, the TS quoted above for the bb channel is larger 
than that found for any of the masses scanned in the t + t~ channel. Our analysis agrees well 
with a standard global LAT ScienceTools analysis using the MINUIT subroutine MINOS [77]. 

The sensitivity of the combined analysis depends most heavily on the observations of the 
dwarf galaxies with the largest J-factors. Specifically, the results of the combined analysis are 
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dominated by Coma Berenices, Draco, Segue 1, Ursa Major II, Ursa Minor, and Willman 1. 
Because the combined analysis is dominated by a few dwarf galaxies with high J-factors, 
any statistical fluctuation in the combined results will necessarily be associated to an excess 
coincident with these dwarf galaxies. We find that the primary contributors to the deviation 
of the combined analysis from the null hypothesis are the ultra-faint dwarf galaxies Segue 1, 
Ursa Major II, and Willman 1. We find a comparable y-ray excess coincident with Hercules 
and slightly smaller excesses coincident with Sculptor and Canes Venatici II; however, the 
relatively small J-factors of these dwarf galaxies limit their contribution to the combined 
analysis. On the other hand, the lack of any y-ray excess coincident with Coma Berenices, 
Draco, or Ursa Minor decreases the significance of the deviation in the combined analysis. 
Thus, we find no significant correlation between the J-factor and the y-ray TS of individual 
dwarf galaxies. We investigate the impact of removing the ultra-faint dwarf galaxies from 
the combined analysis at the end of Section VI. 8 

VI. SYSTEMATIC STUDIES 

In this section we briefly summarize the process for determining the impact of systematic 
uncertainties in the LAT instrument performance, the diffuse interstellar y-ray emission 
model, and the assumed dark matter density profile of the dwarf spheroidal galaxies. A 
quantitative summary of these results can be found in Table VIII. Additionally, we discuss 
the analysis of random blank-sky locations to validate the sensitivity and significance as 
quoted in Section V. Finally, we consider the impact of removing three ultra- faint dwarf 
galaxies from the combined analysis. 

We follow the “bracketing IRF” prescription of Ackermann et al. [28] to quantify the 
impact of uncertainties on the instrument performance. We re-analyze the LAT data us- 
ing two sets of IRFs: one meant to maximize the sensitivity of the LAT to weak y-ray 
sources, and the other meant to minimize it, consistent with the uncertainties of the IRFs. 
When creating these bracketing IRFs, we utilize constraints on the systematic uncertainty 
from Ackermann et al. [28]: 10% uncertainty on the effective area, 15% uncertainty on the 
PSF, and a 5% uncertainty on the energy dispersion. The maximized IRFs increase the 
effective area, decrease the PSF width, and disregard the instrumental energy dispersion. 

8 The results discussed in Section V, including bin-by-bin likelihood functions, are available in a machine- 
readable format at: http : //www-glast . Stanford. edu/pub_data/713/. 
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Conversely, the minimized IRFs decrease the effective area, increase the PSF width, and 
include the instrumental energy dispersion. For each IRF set we re-run the full LAT data 
analysis to fit the normalizations of all background sources. The resulting impact on the 
combined limits for the dark matter annihilation cross section are shown in Table VIII. 
Additionally, it should be noted that discrepancies between the IRFs (derived from Monte 
Carlo simulations) and the true instrument performance could introduce spectral features at 
the ~2-3% level. Because the putative dark matter spectra are more sharply peaked than 
the diffuse background components, spectral features induced by the IRFs could masquer- 
ade as small positive signals. The impact of these effects is not large enough to introduce 
a globally significant signal, but may be partially responsible for the observed discrepancy 
between simulations and LAT data. 

To quantify the impact of uncertainties on the interstellar emission modeling, we create 
a set of 8 alternative diffuse models [78]. Following the prescription of Ackermann et al. 
[79], we generate templates for the Hi, CO, and inverse Compton (IC) emission using the 
GALPROP cosmic-ray (CR) propagation and interaction code. 9 To create 8 reasonably extreme 
diffuse emission models, we vary the three most influential GALPROP input parameters as 
determined by Ackermann et al. [79]: the Hi spin temperature, the CR source distribution, 
and the CR propagation halo height. We simultaneously fit the spectral normalizations of 
the GALPROP intensity maps, an isotropic component, and all 2FGL sources to 2 years of 
all-sky LAT data. Additionally, a log-parabola spectral correction is fit to each model of the 
diffuse y-ray intensity in various Galactocentric annuli. When performing the fit, we include 
geometric templates for Loop I [80] and the Fermi Bubbles [71]. The template for Loop I 
is based on the model of Wollcben [81], while the template for the Fermi Bubbles is taken 
from de Palma et al. [78]. We re-analyze the LAT data surrounding the dwarf spheroidal 
galaxies with each of the 8 alternative models, and we record the most extreme differences in 
the upper limits in Table VIII. Because these diffuse models are constructed using a different 
methodology than the standard LAT interstellar emission model, the resulting uncertainty 
should be considered as an estimate of the systematic uncertainty due to the interstellar 
emission modeling process. While these 8 models are chosen to be reasonably extreme, 
we note that they do not span the full systematic uncertainty involved in modeling the 
interstellar emission. We also note that these 8 models do not bracket the official LAT 

9 http : //galprop . stanford.edu 
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interstellar emission model since the methodologies used to create the models differ. We 
find that variations between the 8 alternative diffuse models have a small impact on the 
combined limits (Table VIII); however, they can increase the TS of the low-mass excess by 
as much as 50%. 

While stellar kinematic data currently do little to constrain the inner slope of the dark 
matter density profile in dwarf spheroidal galaxies, the modeling in Section IV B confirms 
that the integrated J-faetor within an angular radius of 0.5° is relatively insensitive to the 
shape of the inner profile [68]. We find that performing the combined analysis assuming a 
Burkert profile for the dark matter distribution in dwarf galaxies increases the cross section 
limits by ~ 15%. Additionally, we examine the impact of changing the spatial extension 
of the NFW profiles used to model the dwarf galaxies. Taking the ±1 a value of the scale 
radius (as determined in Section IV B) leads to a change of < 20% in the LAT sensitivity. 

We utilize random blank-sky locations as a control sample for the analysis of dwarf 
spheroidal galaxies. (A similar procedure was used by Mazziotta et al. [23] to analyze the 
Milky Way dark matter halo.) We choose blank-sky locations randomly at high Galactic 
latitude (|6| > 30°) and far from any 2FGL catalog sources (> 1° from point-like sources and 

> 5° from spatially extended sources). We select sets of 25 blank-sky regions separated by 

> 7° to correspond to the 25 dwarf galaxies. For each blank-sky location, data are selected 
and binned according to Section II, diffuse background sources and the local 2FGL point-like 
sources are fit, and the likelihood analysis is performed. We map each set of random sky 
locations one-to-one to the dwarf galaxies, and we randomize the nominal J-faetors according 
to the uncertainty derived from kinematic measurements. We form a joint likelihood function 
from sets of 15 random locations to conform to the 15 dwarf galaxies used in the combined 
analysis of Section VB. 

In Figure 4 we plot the TS values derived by fitting for a source with a 25 GeV bb spectrum 
at each of the 7500 individual locations, and in Figure 5 we show the expected sensitivity from 
300 combined analyses performed on randomly selected locations. While the TS distribution 
derived from simulations agrees well with asymptotic theorems, it is clear that the TS 
distribution from random blank fields deviates from this expectation. The LAT data are 
known to contain unresolved point-like astrophysical sources which are not included in the 
background model. Additionally, imperfect modeling of the diffuse background emission 
and percent-level inconsistencies in the determination of the instrument performance may 
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contribute to deviations in the measured TS distribution. Each of these systematic effects 
is a plausible contributor to deviations from the asymptotic expectations, suggesting that 
the significance of the dark matter signal hypothesis cannot be naively computed from 
the TS using asymptotic theorems. While the study of random sky locations suffers from 
limited statistics, increasing the number of randomly selected locations would reduce the 
independence of each trial. 

The combined limits presented in Section V include the ultra- faint dwarf galaxies Segue 1, 
Ursa Major II, and Willman 1. While Bayesian hierarchical modeling sets rather tight con- 
straints on these dwarf galaxies as members of the dwarf galaxy population, the stellar kine- 
matic data yield larger uncertainties on the J-faetors of these dwarf galaxies when analyzed 
individually. For example the velocity distribution of Willman 1 appears to be bimodal, and 
it is not yet clear whether the stars are dynamically bound to this object [55]. Similarly, the 
photometry of Ursa Major II shows that this object is elongated; however, the major axis of 
Ursa Major II does not point in the direction of the Galactic center and it is unclear if this 
elongation should be interpreted as a sign of tidal disruption [82], Removing these three 
ultra-faint dwarf galaxies has a ±20% impact on the combined limits for soft annihilation 
spectra (e.g., bb models or t + t~ models with tttdm ~ 100 GeV) and weakens the limits by 
a factor of ~2 for hard annihilation spectra (e.g., t + t~ models with tudm > 500 GeV); 
however, it does not significantly impact the constraint on the thermal relic cross section 
at low dark matter masses (Figure 6). Moreover, as mentioned in Section V, these ultra- 
faint dwarf galaxies contribute heavily to the TS excess observed in the combined analysis. 
The removal of all three ultra- faint galaxies reduces the maximum excess from TS = 8.7 to 
TS < 0.5. Conversely, forming a composite analysis from just the three ultra-faint dwarf 
galaxies yields TS ~ 10, the significance of which must be decreased by an additional trials 
factor. 


VII. BAYESIAN ANALYSIS 

As a further check of the maximum likelihood analysis, we perform a complementary 
Bayesian analysis based on the technique of Mazziotta et al. [23]. This approach derives 
the diffuse background empirically from annuli surrounding the dwarf galaxies and incor- 
porates the instrument response through an iterative Bayesian unfolding [31, 83, 84], The 
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Bayesian method discussed here improves upon that of Mazziotta et al. [23] by folding the 
dark matter signal spectrum with the IRFs and incorporating information from all energy 
bins when reconstructing a posterior probability distribution for the dark matter cross sec- 
tion. We apply this method to our 4-year data sample and derive upper limits on the dark 
matter annihilation cross section from a combined analysis of the 15 dwarf galaxies selected 
in Section VB. We find that the results of the Bayesian and maximum likelihood analyses 
are comparable despite intrinsically different assumptions and methodologies. 

We begin by selecting signal and background regions centered on each of the dwarf galax- 
ies. Signal regions are defined as cones of 0.5° angular radius, 10 while background regions are 
defined as annuli with inner radii of 5° and outer radii of 6°. When defining the background 
regions, all data within 3° of 2FGL catalog sources are masked. Data are binned in energy 
according to the prescription in Section III. 

Within each energy bin j we denote the number of counts in the signal region as rij and 
the number of counts in the background region as rrij. We assume that the probabilities 
of measuring rij counts in the signal region and rn 3 counts in the background region are 
both Poisson distributed with expected values Sj + cbj and bj , respectively. Here, s 3 is 
the expected number of signal counts (in the signal region), bj is the expected number of 
background counts (in the background region), while c is the ratio of solid angles between 
the signal and background regions. 11 The values of Sj are tied across the energy bins by the 
shape of the putative dark matter spectrum and folded with the LAT instrument response. 
For each dwarf, the response of the LAT instrument is evaluated from full Monte Carlo 
simulations of the detector, which take into account the energy dispersion. 

Starting from the measured counts in the signal and background regions, we evaluate the 
Bayesian posterior probability density function (PDF), 

p((av), J, b ) oc C((av), J,b\n, m) 

x V{{<rv))V{J)V{b) . 

Here we indicate with n = {rij} and m = {rrij} the counts in the signal and background 
regions and with b = {bj} the expected counts in the background regions. The likeli- 
hood C((av ), J, b | n, m) is the product of Poisson probabilities across all energy bins, while 

10 The signal region is chosen to correspond to the solid angle used for J-factor calculations (Section IV B). 

11 Point-source masking is accounted for when calculating the ratio of solid angles between the signal and 
background regions [23]. 
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V({crv}), V(J), and V(b) are the prior PDFs for the cross section, J-factor, and background 
counts, respectively. Prior PDFs for the cross section and background counts are assumed to 
be uniform, while the prior PDF for the J-factor is assumed to be log-normal with mean and 
variance taken from Table I. The posterior PDF on ( av } is evaluated from Equation (15), 
marginalizing the joint PDF over the nuisance parameters b and J . Upper limits on (av) 
are calculated by integrating the posterior PDF. We perform a combined analysis of the 15 
dwarf galaxies in Section 111 using a similar technique, weighting the data associated with 
each source by its J-factor [Equation (15)]. 

We validate the performance of the Bayesian analysis by applying it to random blank 
holds in the high-latitude data chosen in a manner similar to that described in Section VI. 
Similar to the maximum likelihood analysis, we find that the observed cross section limits 
for low dark matter masses are slightly higher than the median expected from random holds, 
but are consistent within statistical fluctuations. The Bayesian analysis presented here is 
signihcantly more sensitive to high-mass dark matter models than that presented in Mazz- 
iotta et al. [23] due to the utilization of all energy bins when testing putative dark matter 
signal spectra. However, the upper limits derived from the Bayesian analysis are slightly 
less constraining than those derived from maximum likelihood analysis. The differing results 
may be ascribed to alternative models of the signal (i.e., finite aperture vs. spatial prohle) 
and of the background (i.e., local annulus vs. global template). Additionally, differences 
between the delta-log-likelihood method and the Bayesian method for setting upper limits 
will lead to more conservative Bayesian limits in the low-counts regime. Despite differ- 
ences in background modeling, the treatment of the LAT instrument performance, and the 
methodology for setting upper limits, the Bayesian and maximum likelihood analyses yield 
comparable results (Figure 7). 


VIII. CONCLUSIONS 

We have reported on 4-year y-ray observations of 25 dwarf spheroidal satellite galaxies 
of the Milky Way. No significant y-ray excess was found coincident with any of the dwarf 
galaxies for any of the spectral models tested. We performed a combined analysis of 15 dwarf 
galaxies under the assumption that the characteristics of the dark matter particle are shared 
between the dwarfs. Again, no globally significant excess was found for any of the spectral 


23 



models tested. We set 95% CL limits on the thermally averaged dark matter annihilation 
cross section incorporating statistical uncertainties in the J-factors derived from fits to stellar 
kinematic data. These limits constrain the dark matter annihilation cross section to be less 
than 3 x 10 -26 cm 3 s -1 for dark matter particles with a mass less than 10 GeV for the bb 
channel and a mass less than 15 GeV for the uu and r + r~ channels. Our limits on the 
dark matter annihilation cross section extend to dark matter masses of 10 TeV, although 
Imaging Air Cherenkov Telescopes present more constraining limits at high dark matter 
masses (Figure 8) [85-87]. 

The analysis presented here greatly improves our understanding of the statistical and 
systematic issues involved in the search for dwarf spheroidal galaxies. Comparing our results 
to the 2-year maximum likelihood study of 10 dwarf galaxies by Ackermann et al. [21], 
we find that the current limits are a factor of ~ 2 weaker for soft spectral models and a 
factor of > 2 stronger for hard spectral models. These changes can be attributed to updated 
J-factors, the increased photon energy range, the inclusion of more dwarf spheroidal galaxies 
(specifically Willman 1), and a reclassification of events when moving from Pass 6 to P7REP. 
As expected, we find that 4 years of P7REP data yield more constraining limits than 2 years 
of P7REP data. 

The largest deviation from the null hypothesis occurs for soft y-ray spectra and is fit by 
dark matter in the mass range from 10 GeV to 25 GeV annihilating to bb and has TS = 8.7. 
Control studies of random blank sky locations suggest that the global significance of this 
excess is p ^ 0.08. This deviation can also be fit by lower mass dark matter (~5GeV) 
annihilating through harder channels (i.e. , t + t~). Much interest in the low- mass WIMP 
regime has been triggered by y-ray studies of the Galactic center [7, 9] and by recent direct 
detection results [88]. While suggestive, the analysis of dwarf spheroidal galaxies presented 
here remains consistent with the null hypothesis within statistical fluctuation. Constraints 
on the dark matter annihilation cross section derived here are found to be robust against 
the systematic uncertainties considered; however, systematics associated with the diffuse 
modeling are found to more significantly impact the TS derived for the low-mass excess. 

Upcoming improvements to the LAT instrument performance and continued data taking 
will lead to increased sensitivity for dark matter annihilation in dwarf spheroidal galaxies. 
The novel event reconstruction and event selection of Pass 8 will be especially important 
as it will both increase the LAT sensitivity to point-like sources and help mitigate sys- 
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tematic effects present in P7REP [89]. Another exciting prospect is the discovery of new 
dwarf spheroidal galaxies. Over the past ten years, data from the Sloan Digital Sky Sur- 
vey (SDSS) [90] have been used to roughly double the number of known dwarf spheroidal 
galaxies, while only covering ~25% of the sky [91]. Thus, the next generations of deep, 
wide-held photometric surveys (e.g., PanSTARRS [92], Southern Sky Survey [93], DES [94], 
and LSST [95]) are expected to greatly increase the number of known Milky Way dwarf 
spheroidal satellite galaxies [96] . The discovery and characterization of new dwarf spheroidal 
galaxies could greatly improve the LAT sensitivity to dark matter annihilation in this class 
of objects. 
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X. FIGURES & TABLES 



FIG. 1. Known dwarf spheroidal satellite galaxies of the Milky Way overlaid on a Hammer-Aitoff 
projection of a 4- year LAT counts map (E > 1 GeV). The 15 dwarf galaxies included in the 
combined analysis are shown as filled circles, while additional dwarf galaxies are shown as open 
circles. 
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Energy (MeV) 


FIG. 2. Histogram of the bin- by-bin LAT likelihood function used to test for a putative 7 - 
ray source at the position of the Draco dwarf spheroidal galaxy. The bin-by-bin likelihood is 
calculated by scanning the integrated energy flux of the putative source within each energy bin 
(equivalent to scanning in the spectral normalization of the source). When performing this scan, 
the flux normalizations of the background sources are fixed to their optimal values as derived from 
a maximum likelihood fit over the full energy range. Within each bin, the color scale denotes the 
variation of the logarithm of the likelihood with respect to the best fit value of the putative source 
flux (truncated at — Alog£ > 10). Upper limits on the integrated energy flux are set at 95% CL 
within each bin using the delta-log-likelihood technique and are largely independent of the putative 
source spectrum. 
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TS 

FIG. 4. Distribution of TS values from individual fits of a 25 GeV bb annihilation spectrum to the 
null hypothesis generated from 50000 realistic Monte Carlo simulations and 7500 random blank-sky 
locations at high latitude in the LAT data. The distribution of TS values derived from simulations 
of individual ROIs is well matched to the expectations from the asymptotic theorem of Chernoff 
[70], while the distribution derived from random sky positions shows an excess at large TS values. 
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FIG. 5. Constraints on the dark matter annihilation cross section at 95% CL derived from a 
combined analysis of 15 dwarf spheroidal galaxies assuming an NFW dark matter distribution 
(solid line). In each panel bands represent the expected sensitivity as calculated by repeating 
the combined analysis on 300 randomly-selected sets of blank fields at high Galactic latitudes in 
the LAT data. The dashed line shows the median expected sensitivity while the bands represent 
the 68% and 95% quantiles. For each set of random locations, nominal J-factors are randomized 
in accord with their measurement uncertainties. Thus, the positions and widths of the expected 
sensitivity bands reflect the range of statistical fluctuations expected both from the LAT data and 
from the stellar kinematics of the dwarf galaxies. The most significant excess in the observed limits 
occurs for the bb channel between lOGeV and 25GeV with TS = 8.7 (global p-value of p ~ 0.08). 
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"Idm (GeV) 

FIG. 6. Constraints on the dark matter annihilation cross section (t + t~ channel) at 95% CL 
derived from a combined analysis excluding three ultra- faint dwarf galaxies: Segue 1, Ursa Major II, 
and Willman 1 (solid line). The expected sensitivity is also calculated excluding these three ultra- 
faint dwarf galaxies and is represented in the same manner as in Figure 5. 
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FIG. 7. Comparison of constraints on the dark matter annihilation cross section (t + t~ channel) 
derived from the combined maximum likelihood and the combined Bayesian analyses of 15 dwarf 
spheroidal galaxies. The expected sensitivity for the maximum likelihood analysis is represented 
similarly to Figure 5. The observed Bayesian limits are consistent with the expected Bayesian 
sensitivity bands (not shown), which are likewise higher than those of the maximum likelihood 
analysis. 
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FIG. 8. Comparison of constraints on the dark matter annihilation cross section (bb channel) 
derived from the LAT combined analysis of 15 dwarf galaxies (assuming an NFW profile), 48-hour 
observations of Segue 1 by VERITAS (assuming an Einasto profile) [85], and 112-hour observations 
of the Galactic center by H.E.S.S. (assuming an Einasto profile) [86]. In the interest of a direct 
comparison, we also show the LAT constraints derived for Segue 1 alone assuming an Einasto 
dark matter profile consistent with that used by VERITAS [85]. For this rescaling, the J-factor 
of Segue 1 is calculated over the LAT solid angle of AH ~ 2.4 x 10~ 4 sr and yields a rescaled 
value of 1.7 x 10 19 GeV 2 cm -5 sr (uncertainties on the J-factor are neglected for comparison with 
VERITAS). 
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J-factors are calculated over a solid angle of Afl ~ 2.4 x 10 4 sr. 
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Willman 1 7.27e-25 8.67e-25 2.11e-24 4.74e-24 1.18e-23 4.44e-23 1.24e-22 3.74e-22 

Combined 1.56e-25 1.84e-25 4.12e-25 8.37e-25 1.75e-24 4.52e-24 1.10e-23 3.27e-23 



TABLE VIII. Systematic uncertainties on the maximum likelihood analysis decomposed into con- 
tributions from the IRFs, the diffuse modeling, and the spatial extension of the dark matter profile. 
Entries represent the percentage change to the combined upper limits on the dark matter annihi- 
lation cross section from varying each component individually. 




lOGeV 

100 GeV 

1000 GeV 

10000 GeV 


IRFs 

+14%/-12% 

+12%/-10% 

+ll%/-9% 

+ll%/-9% 

e + e“ 

Diffuse 

+3%/-4% 

+3%/-3% 

+1%/-1% 

+1%/-1% 


Extension 

+7%/-5% 

+17%/-11% 

+ll%/-6% 

+10%/-6% 


IRFs 

+15%/-12% 

+12%/-10% 

+ll%/-9% 

+ll%/-9% 


Diffuse 

+4%/-5% 

+3%/-4% 

+2%/-l% 

+1%/-1% 


Extension 

+7%/-5% 

+15%/-10% 

+ll%/-6% 

+10%/-6% 


IRFs 

+15%/-13% 

+12%/-10% 

+13%/-11% 

+12%/-10% 

T + T - 

Diffuse 

+5%/-5% 

+l%/-5% 

+1%/-1% 

+0%/-l% 


Extension 

+6%/-4% 

+14%/-9% 

+13%/-7% 

+6%/-3% 


IRFs 

+15%/-14% 

+14%/-12% 

+12%/-10% 

+ll%/-9% 

uu 

Diffuse 

+9%/-4% 

+3%/-4% 

+3%/-4% 

+2%/-3% 


Extension 

+4%/-3% 

+9%/-7% 

+12%/-8% 

+12%/- 7% 


IRFs 

+15%/-13% 

+14%/-12% 

+12%/-10% 

+ll%/-9% 

bb 

Diffuse 

+10%/-3% 

+4%/-5% 

+2%/-4% 

+2%/-3% 


Extension 

+3%/-2% 

+8%/-6% 

+12%/-8% 

+12%/- 7% 


IRFs 

- 

+14%/-12% 

+12%/-10% 

+13%/-11% 

W + W~ 

Diffuse 

- 

+4%/-4% 

+2%/-4% 

+1%/-1% 


Extension 

- 

+8%/-6% 

+13%/-8% 

+14%/-8% 
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